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Abstract 

The orthogonal polynomials with recurrence relation 

(An + Hn ~ z) F n (z) = /i„ + i F n+ i(z) + A n _i F n -i(z) 

with two kinds of cubic transition rates X n and /j, n , corresponding to indeterminate 
Stieltjes moment problems, are analyzed. We derive generating functions for these 
two classes of polynomials, which enable us to compute their Nevanlinna matrices. 
We discuss the asymptotics of the Nevanlinna matrices in the complex plane. 



1 Introduction 



The field of indeterminate moment problems applied to birth and death processes has 
been quite active in the past ten years and many explicit examples have been worked out, 
see [3] and the many references therein. Restricting ourselves to the case of polynomial 
transition rates A n and /i n the results obtained dealt mainly with quartic rates [3], [6]. It 
is the aim of this article to show that the same underlying ideas that led successfully to 
the computation of the Nevanlinna matrices for the quartic rates can be applied to some 
cubic rates, leading to some explicit integral representations for their Nevanlinna matrix 
elements. 

The plan of the article is the following. In Section 2 we will recall some basic relations 
between birth and death processes and orthogonal polynomial theory. In this article we 
will be concerned with the two processes: 

(PI): A n = (3n + 3c + l) 2 (3n + 3c + 2), /i n = (3n + 3c - l)(3n + 3c) 2 (l - 5 n0 ), 

and 

(P2) : A n = (3n + 3c+ l)(3n + 3c + 2) 2 , fi n = (3n + 3c) 2 (3n + 3c + 1)(1 - S n0 ), 

under the assumption c > 0. In Section 3 some background material useful in the sequel 
is gathered. In Section 4 and 5 we obtain some generating functions which will allow, in 
Section 6 to compute the Nevanlinna matrices for both processes. In Section 7 we analyze 
the asymptotics, in the complex plane, of the Nevanlinna matrix elements. 

2 Birth and death processes versus orthogonal poly- 
nomials 

Birth and death processes are special stationary Markov processes whose state space is N, 
representing for instance some population. We are interested in the time evolution of such 
a population, described by the transition probabilities V m ,n{t) yielding the probability that 
the population goes from the state m at time t — to the state n at time t > 0. This 
evolution is supposed to be governed by 

Vn,n+l(t) = A n t + 0(f), 

V n , n (t) = l-(\ n + ll n )t + 0(t), t->0. 
V n , n -i(t) = /J, n t + 0(t), 

For applications the most important problem is to find V m ,n(t) for given rates A n and fj, n , 
with suitable extra constraints to be described later on. 

From the previous setting one can prove that the transition probabilities have to be a 
solution of the forward Kolmogorov equations 

~~^Pm,n A ra _i 'Prn,n—1 f-^n+l Pm,n+1 (An ~\~ fin) T^m,n- (1) 

The V m , n {t) are assumed to be continuous for small time scales with 

Jim V m ,n{t) = S m , n . (2) 



A representation theorem for 7 3 m , n (t) was proved by Karlin and MacGregor in [5] which 
links birth and death processes and orthogonal polynomials theory Let us define the 
polynomials F n (x) by the three-terms recurrence relation 

(\ n + /j, n - x)F n (x) = n n+1 F n+1 {x) + A„_iF„_!(x), n > 1, (3) 

with the initial conditions 

F (x) = 1, F l{ x) = Wo -" 

/ii 

Let us define 

A Ai • • • A n _i 

7T = 1, Tin = , n > 1. 

A*iA*2 • • • 

If the positivity conditions 

A„ > 0, n > 0, and /x = 0, /x n > 0, n > 1 (4) 
are fulfilled, then there is a positive measure tp for which 

Vm,n(t^) = — I e- xt F m (x)F n (x) di>{x). (5) 

Jsupp(ip) 

Then the initial condition (2) is nothing but the orthogonality relation 

F m (x)F n (x) d^{x) = <J m>n . 



^"m J supp(V>) 

Such a measure has well-defined moments 



„= / x n dij(x), n = 0, 1,... 

•/ supp(</>) 



'supp(V') 

If supp(^) C H. this is a Hamburger moment problem and if supp(^) C [0, +oo[ this is 
a Stieltjes moment problem. In the event that the measure ip is not unique we speak of 
indeterminate Hamburger (or indeterminate Stieltjes) moment problems, indet H or indet 
S for short. Stieltjes (see [1]) obtained the necessary and sufficient conditions for a moment 
problem to be indet S 

Kn <oo, ^ < co. (6) 

n>0 n>l An7r ™ 

These conditions imply that it is also indet H. 

Let us consider now the two processes to be analyzed in this article. We will denote the 
first one as the process PI, with rates 

A n = (3n + 3c+l) 2 (3n + 3c + 2), fi n = (3n + 3c - l)(3n + 3c) 2 (l - 5 n0 ), n > 0, (7) 

and the second one as process P2, with rates 

A n = (3n + 3c+ l)(3n + 3c + 2) 2 , fi n = (3n + 3c) 2 (3n + 3c + 1)(1 - 8 n0 ), n > (8) 



Using the notation (a) n = T(n + a)/T(a) we have for PI the large n behaviour 

n n = 2 = 0(n- 4 / 3 ), -±- = 0{n-^) 

\ (c + l) n J ir n fi n 

and for P2 

((c+l)„) 2 (c + 4/3)„ 1 ' - i„ M „ ^ ' 

These asymptotic estimates show that the conditions (6) are satisfied and therefore the 
two processes are indet S and indet H. 

3 Background material 

In order to describe the Nevanlinna matrix we will need a triplet of elementary functions 
defined by 

3n+l 

n>0 V ' 



It is easy to check the relations 

a[ = (To, 2 = 01, Q = -0 2 , 

(10) 



O (O) = 1, 0i(O)=O, 2 (O) = O. 

These functions are called trigonometric functions of order 3, since they are three linearly 
independent solutions of the third order differential equation 

a'/' + a { = 0, 1 = 0,1, 2. 

Their explicit form is 

a (u) = ^(e- u + e ju + e ]U ) = ^ ( e~ u + 2 cos ( ^-u ] e u/2 



0i («) = i(- e -« + je>» + = I i-e~ u + 2 cos (^u - | ) e"/ 2 ) , j = e^ 3 . 



<7 2 (u) = \{e- u - je>« - Je*) = i ( e "« - 2 cos ( + | ] e"/ 2 



We will need also the following functions 



:n) 



s « = jf <r^ = ^1(1^ = ^' (12) 

Observing that 6>(t) is continuous, decreasing and concave for t e [0, 1] gives the bounds 



4 First generating function 



We will consider, for c > 0, the slightly more general rates than the ones defined in (7): 
A n = (3n + 3c+ l) 2 (3n + 3c + 2), n>0 



H n = (3n + 3c — l)(3n + 3c) 2 , 



n > 1 



(14) 



where /xo is taken as a free parameter, not necessarily equal to (3c — l)(3c) 2 . 

We will denote by F n (z;c, Ho) the polynomials with recurrence relation (3) and the 
rates (14). Obviously the polynomials corresponding to process PI are recovered as the 
limiting values F n (z; c, 0). 

To get most conveniently a generating function for them, it is useful to define a triplet 
of functions d 3n+ i((), I — 0, 1, 2 for n — 0, 1, ... by the recurrence relation 



n > 



with the initial values 



d 3n +l = _ C d 3n + ^3n-2, 
d?,n+2 = —(d^n+l, 

K d 3n +s = —( d 3n+2 + A n d 3n , 



(a) 
(b) 
(c) 



(15) 



d. 



1 



do 



and the definition ( = z 1 ^. Note that d 3n (() are polynomials with respect to z. 
Let us begin with 

Proposition 1 The polynomials F n and d n are related by 

(c+l/3) n d 3 n(C) 



F n (z;c,no) = (3c)! 



(16) 



(c+l) n (3n + 3c)!' 
where we use the notation (a)\ = F(a + 1) for a > 0. 
Proof : 

Let us define M n (z) = d 3n ((). Using (15c) and (15b) in (15) we have for n > 1 

M n+ i = d 3n+3 = -( d 3n+2 + A n d 3n = A n M n + C 2 <^3n+i = (A„ - z)M n + ( 2 (( d 3n + d 3n+1 ). 

Then we use (15a) and (15b) with the shift n — > n — 1 to get 

M n+1 = (A n - z)M n + fj^ C 2 d 3 n-2 = (A n - z)M n - fl n ( d 3n ^. 

The term involving d 3n -\ is disposed of using the (15c) with the shift n — > n — 1. One is 
left with 



M n+ i = (A„-2)M n +/i n (rf 3n -A„_irf3 n _3) = (A n +/j n -20M„-A„_i/j n M n _i, 
The boundary conditions are to be computed separately and give 

M Q = 1, M 1 =d 3 = \o + j2o-z. 



n > 1. 



It is then easy to check that the polynomials F n are related to the M n by 

j? , \ M n{z) d 3n (C) , , 

F n {z; c, no) = = , n>0. (17) 

Using Gauss multiplication formula we have 

(3r l + ,| C)! = 3 3 "(c + l/3) B (c + 2/3) B (c + 1) B = ^ff" ^ • • • n n , 
(3c) I (c+l) n 

and this leads to the desired relation (16). □ 

In view of Proposition 1 we need generating functions for d 3n+ i which we define for the 
variable t - in a suitable neighbourhood of the origin - as 

j.3n+3e+2 

G,(C,t) = E^,(C) c+ , ( = 0,1,2. (18) 

n>0 y ' 

Routine computations, using relations (15) give for these generating functions the linear 
differential system 

y-3c-l 

(l-t 3 )D t G -t 2 G + (G 2 = 



(1 - t 3 ) D t Gi - 2t 2 d + (G 



(3c- 1)!' 
A*o t 3c (19) 



C 2 (3c)!' 
D t G 2 + (G 1 = 0. 

All the factorials involved are well defined in view of the hypothesis c > 0. 
Switching from the functions Gi to new functions Gi defined by 

Go G\ ^- 

G ° = (1 _f3)l/3' Gl = (1 _ £3)2/3' ° 2 = G2 ' ( 20 ) 

the differential system takes the more symmetric form 
(l-t 3 ) 2 / 3 D t G + (G 2 ' 



(3c- 1)!' 



(i-t 3 r/ 3 D t d 1+ cG = ^^(i-^ 3 ) 173 ' (21) 

(\-t 3 yi 3 D t G 2 + (G l = 0. 

Using the variable 9{t) defined in (12) we observe that (1 — t 3 ) 2 ^ 3 D t = D g so that (21) 
becomes an inhomogeneous differential system with constant coefficients: 

DgG + CG 2 = a{6), D e G 1 + (G = ^, D e G 2 + (G 1 = 0. (22) 

This is easily solved for G ; one gets 

G (9) = f a (((9 - v)) a(v) dv + f a ^ 9 ~ v ^ b ( v ) dv. (23) 
Jo Jo s 



From this result we recover Gq(z, t) by coming back to the original variable t and after the 
change of variable v = 9~ l (u) in the integral. Using the notation 0(t, u) — 9(t) — 0(u), we 
conclude to: 



G {z,t)= / a ((e(t,u))- _(l_«S)-2/3 du 

Jo (3c- 1)! 



-Mo 



(24) 



*2(m,u)) U (1 _ u s r i/3 dli . 



C 2 (3c)! 

Taking into account relations (20), (17) and (16) we have on the one hand 

g (z, t) = (i - t 3 y/ 3 g (z, t) = j^J2 rrrk F ^ c ' ^ t3n+3c Q ~ t3 ) 1/3 ' 

^ >' n>0 ^ C l l 6 )n 

and on the other hand G (z,t) given by (24). Gathering all these pieces we end up with 
Proposition 2 The polynomials F n (z; c, /xo) have the generating function 

Yl JTTTP^r Fn ^ c > t3 " +3c ( x ~ t3 ) V3 = 

3c f a ((Q(t,u))u 3c -\l-u 3 r 2 / 3 du + fi f ' U)) u 3c (l - u 3 )^ 3 du. 

Jo Jo C 

This is not quite enough to compute the Nevanlinna matrix; in fact we need 

F(z;c,hq) = ^2 F n(z;c,fj, Q ). (26) 

n>0 

Using the notation 

i»r(/3) 
5(a ^ )= r(^)' 

we will now state: 

Proposition 3 We have the relations 



+„„ /'^(i-Q-v 1 -^)) 

Jo 2 



(27) 



va/id /or c > and 



^ 3 (3c- l)(3c) 



1 c,/i ) 



-/ u 3c - 3 {l-u 3 y l/3 o- {(6{u))du 
Jo 



Mo 5^-2/3,2/3) /i 

fro-(3c-l)(3c) 2 ] 3 _ /3 (C £ (U)) d 

/x £?(c+l/3,2/3) 7 1 J 

(28) 

va/id /or c > 1 . 



Proof : In (25) we set t = r 1//3 , then we multiply both sides by r~ 2 / 3 (l — r)~ 2 / 3 and 
integrate from r = to r = 1. The left hand side integral involves a Eulerian integral and 

we get 

r(c+l/3)r(2/3) ^ , . »,„., t 1 G (z,t) 



r (c + i) E/»fe^o) = 3(30!/'^^ 



The right hand side is a double integral, which, upon interchange of the order of the 
integrations and use of relations (10), gives (27). In this last result, the integral with no a 
function, when expressed in terms of Euler gamma functions simplifies to fio/z. The first 
integral in (27), using relations (10), can be integrated by parts twice; then elementary 
algebra results in (28). □ 

The results obtained so far are sufficient to compute the functions C\{z) and Di(z) in 
the Nevanlinna matrix of process PI. However, to get the full matrix we need also the 
generating function for the dual process in the sense of Karlin and McGregor (KMG for 
short). For the reader's convenience let us recall its definition. 

The correspondence from a process V to its KMG dual V is as follows 

P — {^n, A*n} —> V = {\ n = fi n+ l, fi n = A n }. 

It follows that the dual process of PI will have 

A n = (3n + 3c + 2)(3n + 3c + 3) 2 , Jl n = (3ra + 3c + l) 2 (3n + 3c + 2), n > 0, (29) 

which correspond to the process P2 up to the shift c — > c + 1/3. So we will now work out 
a generating function for the process P2. 

5 Second generating function 

Here again we will consider, for c > 0, the slightly more general rates than the ones defined 
in (8): 

A n = (3n + 3c+ l)(3n + 3c + 2) 2 , n>0 

(30) 

H n = (3n + 3c) (3n + 3c + 1), n > 1 /io > 0, 

where /io is taken as a free parameter, not necessarily equal to (3c + l)(3c) 2 . 

We will denote by G n (z; c, /io) the polynomials with recurrence relation (3) and the rates 
(30). Obviously the polynomials corresponding to process P2 are recovered as G n (z; c, 0). 

In order to avoid repetitions, we will give only the main steps. It is again useful to 
define a triplet of polynomials e 3n+ ;(C), / = 0, 1, 2 by the recurrence relation 



C = z 1/3 n > 



e 3n+l — —( &3n + G3n-2, 

G3n+2 = — C e 3n+l) (31) 
, e 3n+3 — —( e 3n+2 + &3n, 



with the boundary values 



e_i = --, e = 1, =>• e 3 = A + /io - -2- 

By an argument which follows closely the one given in the proof of Proposition 1, we get: 



Proposition 4 The polynomials G n and e n are related by 



n ( \ ,a , iM (c + 2/3) n e 3n (C) 

G n (Z\ C, Un) = OC + 1 ! — r — . 

y ^ ! V ; (c+l)„ (3n + 3c+l)! 



(32) 



We then define the generating functions 

^3n+3c+l 

e 3n+l{Q) ^ 

n>0 

for which we get the differential system 

t 3c-l 



4.6U+6C+1 

g,(C,«)=E^(C) (3n + 3c + or ' = 0,1,2 (33) 



(l-^)A^o-2t 2 // + C^2 = 



(3c -1)! 

D t H x + CHo = (34) 

(1-t 3 ) D t E 2 -t 2 E 2 + C,E x = -^-J—- 
V ; S C (3c +1)! 

Switching from the functions ifj to new functions ifj defined by 

H ° = (i _t°)2 /3 ' Gl = ^' ^ 2 = (1 -t 3 ) 1/3 ' 

and using the variable 6(t) defined in (12) the previous system becomes an inhomogeneous 
differential system with constant coefficients, easy to solve. Combining all this we get: 

Proposition 5 The polynomials G n (z;c, // ) have the generating function 
f^(c + 2/3) n y c 

r V 2 (ce(t,,)) M 3 C+1(i _ M3) . 2/ 3^ 

Jo ( 

(36) 

In fact we need the generating function 

Q(z\ c, no) = E G n(^ c > A*o)- (37) 

n>0 

Let us prove: 

Proposition 6 We have the relation 



5 <^"°> = B ( c + 2/3.1/3) { 3C(3C+ / ^f^^ 1 (1 - «')-"*. 

Jo 2 



(38) 



valid for c > and 

1 - i ^' C ' "»> = B(c- 1/3.1/3) 1!££ f ±il / ' "" 2(1 - « 3 )- 2/V »« ? (")> *' 
+ B (e + 2/3,l/3) " 3C+1(1 " du, 

(39) 

va/id /or c > 1/3. 

Proof : In relation (36) we change the variable t to r defined by t = r 1//3 , then multiply 
both sides by t~ 2 / 3 (1 — t)~ 2 / 3 and integrate from r = to r = 1. The left hand-side is 
merely a Eulerian integral, while the right-hand side is a double integral. Interchanging the 
order of integrations, and using relations (10) one gets (38). The integral which does not 
involve a functions can be expressed in terms of Euler Gamma functions. Then elemetary 
algebra yields 39). □ 

Equipped with these results, let us turn ourselves to the determination of the Nevan- 
linna matrix for the processes PI and P2. 

6 The Nevanlinna matrices 

We will write the first Nevanlinna matrix as 

/ Mz) C^z) \ 
MW= • (40) 

\ B 1 (z) D 1 {z) J 

As shown in [6] one gets simpler results by considering the modified Nevanlinna matrix 

=1 I , (41) 



where 



B^z) D^z) 

M*) - M*) - Mz) = Bl (z) - ^M, -1 = £ J-. 

a a a ^—f /i n 7r n 

n>l 

Let us begin with the computation of the modified Nevanlinna matrix for the process 



PI. 

6.1 The Nevanlinna matrix for PI 

We have first 

Proposition 7 The modified Nevanlinna matrix Afi(z) of process PI can be expressed in 
terms of the generating functions T and Q as 

Mz) = T- S{z; c + 1/3, A ) - "r - Q{z- c + 4/3, Ai), d(z) = 1 - — T{z\ c + 1, m), 
Ao Ai /il 

= -l + ^£(z;c+l/3,A ), £>i(z) =2^;c,0). 



Proof : 

We use successively the relations proved in Lemma 6 of [6]. For the reader's convenience 
we will recall these relations. We have for the first element 

D 1 (z)=zJ2 F n (z;c,0) = zF(z;c,0), (42) 

n>0 

upon use of (26). The second element C±(z) is given by 

C 1 {z) = l-^ ) Y J F£\z ] c), (43) 

A'O n>0 

where the polynomials Fn\z; c) have the shifted rates 

Ai 1} = A n+1 , rf'=^i =► FW(z;c) = F n (z;c+l,iJi 1 ). (44) 

Using (26) we conclude to 



C 1 (z) = l-—F(z;c+l,n 1 ). (45) 
Hi 



The third element is given by 

B 1 (z) = -l + ^-Y t F n (z;c), (46) 

^ n>0 

where the F n (z; c) are the KMG duals of F n (z; c), given here by 

F n (z; c) = G n (z; c + 1/3, A ) B^z) = -1 + ^G(z; c + 1/3, A ), (47) 

where the last equality follows from (37). The last element is given by 

Ai^) = ^J2^c)-^J2^ 1) (z;c). (48) 

We can write, using (37) 

F^\z; c) = G n (z; c + 4/3, Ai), Jli — Ai, (49) 

from which the proposition follows. □ 

The Nevanlinna matrix follows quite easily now: 

Proposition 8 The Nevanlinna matrix of the process PI, with rates 

\ n = (3n + 3c + l) 2 (3n + 3c + 2), /i n = (3n + 3c - l)(3n + 3c) 2 (l - 5 n0 ), 

with c > 0, is given by 



B(c + 1,1/3) (3c +1) 7 ' v 7 C 2 

Cl= g( C+ l/3,2/3) f^l-^ ^oW)^, 

A = ' , , 3c, f\^\l-u*)- 2/ * aiiCd{U)) du, 

£?(c+l/3,2/3) 7o C 



(50) 



where ( = z 1 ^. 



Proof : 

The matrix element D\ follows from Proposition 7 and (27). The matrix element C\ 
follows from Proposition 7 and (28). The matrix element B\ follows from Proposition 7 
and (28). To compute A\ we use Proposition 7 and (39) to get first 



zAi 



1 



B(c+ 1,1/3) 3c + 1 J 
which is nothing but 

7 3 1 



Zcu^- 1 - (3c + l)w 3c+2 (1 - u 3 )- 2/3 <T (C6(u)) du. 



B(c + 1,1/3) 3c +1 J 



D, 



u 



3c (l _ u 3 ) l/3 



(Jo((0(u)) du. 



An integration by parts gives the required result in Proposition 8. □ 

Let us observe that the matrix elements Di and Bi can be simplified for c > 1/3. Using 
relations (10) one realizes that an integration by parts is then possible, leaving us with the 
matrix elements 



B, 



B(c+ 1,1/3) 3c +1 J 



u-^ldu, 



(51) 



B(c+ 1/3,2/3) 



u 3c ( a 2 ((9(u)) du. 



The special cases where c = and c = 1/3 lead to considerable simplifications for B and 
D given in [4]. 



6.2 The Nevanlinna matrix for P2 

There is no need to give again the detailed proofs, since everything proceeds as for process 
PI. Beware that now 

A n = (3n + 3c + l)(3n + 3c + 2) 2 , pi n = (3n + 3c) 2 (3n + 3c + 1)(1 - S n0 ). 

We have first 



Proposition 9 The modified Nevanlinna matrix A/^z) of process P2 can be expressed in 
terms of the generating functions T and Q as 

A 2 {z) = 1- J={z- c + 2/3, A ) - ^ T{z- c + 5/3, Ax), C 2 (z) = 1 - — Q(z; c + 1, ^), 
B 2 (z) = -1 + ^- ^(z; c + 2/3, A ), £> 2 (*) = ^ ^(^; c, 0). 

Combining this result with the explicit forms of these generating functions, and upon 
integrations by parts, we get 



Proposition 10 The Nevanlinna matrix of the process P2, with rates 

\ n = (3n + 3c + l)(3n + 3c + 2) 2 , /i n = (3n + 3c) 2 (3n + 3c + 1)(1 - <5 n0 ), 
where c > 0, gwen 6j/ 



3 i r 7j3c ^(W) 

B(c+ 1,2/3) (3c + 2) 7 C 2 



(52) 



^ + 2) 7o 

C 2 = ? " „(3c + l) f\-^pi d u, 
B(c + 2/3,1/3) 1 Wo C 

^ - B (c + 2/3,l/3) + ^ / 1 ""^ " " 3 >- 1/3 ^ *" 

Here too, the c — ► limit gives again simplifications of the matrix elements B and D, see 
?!]• 

We will now examine the asymptotics of the entire functions appearing in the Nevan- 
linna matrix. 

7 Asymptotics of the Nevanlinna matrices 

Three quantities are essential to describe the large \z\ behaviour of an entire function A(z), 
with Taylor series 

A{z)=Y J UA)z n . (53) 

n>0 



The order p A is defined as 



nlnn 

" J = limsup iHU^)H' (54) 



If it is finite, we can define the type a A as 



a A = — hmsup n\UA)\ pA/n , (55) 

the Phragmen-Lindelof indicator h A ((f)) being defined by 

h A {<t>) =limsup ln l A ( re ^)l ; 0g[O,2tt]. (56) 

The knowledge of the indicator gives the type via the relation 

a A = sup h A (4>). (57) 

0G[O,2tt] 



As a preliminary remark, let us observe that all the matrix elements of the Nevanlinna 
matrices have the generic structure 



JV,(z) = f /(«)£,(*, u) du, f(u) = u a (l - u 3 ) b , 1 = 0,1, 2, 
Jo 



(58) 



possibly up to a single factor of z, appearing in D 1 and D 2 - We will not care about this 
factor since it does not change the order, type and Phragmen-Lindelof indicator. 
The entire functions E\ which appear have the structure 

= ^ = £<^-<M*» (59) 

By inspection, we see that the possible values of the parameters are 

a = 3c, 3c -1, 6 = 0,-1/3,-2/3 
so that in any case we have 

a > -1 and b > -2/3, (60) 
and these conditions ensure that / is integrable over [0, 1]. 

Let us begin with 

Proposition 11 The order of the entire functions N t (z), I = 0, 1, 2 is p = 1/3. 
Proof : 

Since E[(z,u) are entire functions of z, uniformly in u G [0, 1], we can integrate term by 
term in relation (59). This gives 

N t (z) = 5>,„* B , 6,n = T^^y f f(u)e 3n+l (u) du. (61) 

n>0 [6n + L).J 

From (61) we have 

hn= T f{u)P n+l {u)du. (62) 
Jo 



n In n n ln n n ln n ' 



In the large n limit, using Stirling formula, this relation yields 

— = 3 — hm sup — . 

p n->oo n In n 

To prove that the second term vanishes in that limit, (because of the logarithm) we need 
an upper and a lower bound for the integral. Using the inequalities given in (13) we get 

dl n+l [ u a (l - u) b+l+Zn du < I hn < 6l n+l M, M= f f(u) du. 
Jo Jo 

It follows that the logarithm of the upper and lower bounds behave, for large n, respectively 

as 

(3n + Z)ln# - (a + l)mn, (3n + I) ln 9 + ln M, 



and these imply 

hi h n 
hmsup — - — — = 0, 

„^oo n In n 

and l/p — 3. □ 

It follows that all the matrix elements of the two Nevanlinna matrices have the same 
order 1/3. Let us state this result as: 

Proposition 12 All the matrix elements of the Nevanlinna matrices given by Proposition 
8 and Proposition 10 have one and the same order p = 1/3. 

Let us now determine the Phragmen-Lindelof indicator. We first need to prove the 
following lemma 

Lemma 1 Under the hypotheses (60) we have, fort — > +oo ; the asymptotic behaviour 

I(t) ee jf 1 u\l - e -«*<«> du = e -*(-+D* + O , z = e» (63) 

provided that <fi e] — 7r/2, +7r/2[. 
Proof : 

The correspondence w — > x = 6(u) is a continuous bijection from [0, 1] — > [0, ^o]- It follows 
that its inverse function 6~ l is also a continuous bijection. Taking i as a new variable, the 
integral I(t) can be written as 

I(t) = xtx) e~ tz * dx, xtx) = [6-\x)ni - {e-\x)f] b+2 '\ (64) 

In view of the hypotheses (60) the parameter 6 + 2/3 is positive and x( x ) * s continuous 
over ]0, ^ ]- Let us split the integral into two pieces: 

X (x) e~ tzx dx, I 2 (t) = ^ X (x) e~ tzx dx, 
Since x is continuous over [e, 9 ], we have 

g— t(Rc z)e g— t(Rez)6o 

\h(t)\ < sup x(x 



xe[e,e ] t(Rez) 

and these terms are vanishing exponentially for large t, so they will be negligible when 
compared to inverse powers of t. 

Next let us consider I\. For sufficiently small e we can write x( x ) = x a + 0(x a+z ). So 
we get 

h(t)= / [x a + 0(x a+3 )]e- tzx dx = / / .... (65) 

J0 J0 Je 

The second integral in (65) is again exponentially small for large t since we have 

/oo rco poo 

x a e~ tzx dx < / x a e- t(Rez)x dx < e" e * (Rc2)/2 / x a e~ t(Rez)x/2 dx, 



and similarly for the term containing 0(x a+ ). The further bounds 



x a e tzx dx 



<e -et{Kcz)/2 I x a e -t(Rez)x/2 dx 



T(a+1) 



3 -ei(Rez)/2 



[t(Rez)/2] a 



+1 



show that this term is also with an exponential decrease. The first integral in (65) is well 
known to give 



Jo 



x a e~ tzx dx 



r(q + l) 



i{a+l)<t> 



So we end up with 1 

+ f 0(x a+3 )e- tzx dx. 
Jo 

For the last piece its modulus is bounded by 

f x a+3 e -t(Rcz)x dx< f x a+3 e -t(Rcz)x dx = q (A 

Jo ~ Jo V a+A J 

so the lemma is proved. □. 

In order to get the indicator let us point out that in the Nevanlinna matrix elements, 
only the entire functions ai((6(u)) / '( l , / = 0, 1, 2 do appear. Since the factor ( l does not 
change the indicator (see the definition (56)), it is sufficient to deal with the o\ and from 
relation (11) these can be written as linear combinations of exponentials. Considering a 
generic matrix element of the form 



N(z)= f f(u)[a e-^ + a. 
Jo 



e K0(u) + a _ e 



z = C 3 , (66) 



we will prove: 

Proposition 13 All the functions N have as Phragmen-Lindeldf indicator 

'(f) — 7T N 



h(4>) = 9 cos 



<pe [0,2tt]. 



(67) 



Proof : 

Let us use the notations 

z = e i(p/3 = c + is , 



z + = e i{<t,+n)/? ' = C+ + is 



z _ _ e i(4>-n)/3 = c _ + - s _ ? 



let us determine the Phragmen-Lindeldf indicator first for e]0,27r[. 
It is elementary to check that the following inequalities hold: 



0e 
0e 



7T 3lT 

2' T 
3tt 



c_ > c + > -c < 

c_ > y/3/2 -c < 0, c + < 



(68) 



-,2tt 



c_ > -c > 



c+ < 



1 All the equalities are understood up to exponentially small terms omitted. 



The asymptotic behaviour of the integral must be analyzed separately for these three cases. 
Let us begin with the simplest case, when e]f , ^[. One has 



N(z) 



J(r,0) 



with 

J(r,0) = a_ f f{u)i 
Jo 



1 ' i e{u) 



du 



+ e~ z - r 1/:% > f f(u)[a + e z + rl/3d ^ + a^ 11 '^} du. 
Jo 



We have merely used the relation 6{u) = 9 — 6{u) in the piece involving a_. Since the real 
parts of z + and — z Q are negative, the absolute value of the last two integrals is bounded. 
Furthermore, since the real part of z_ is positive, for large r the factor e ~ z - rl/3e o h as an 
exponentially decreasing absolute value. So we conclude to 



lim sup 



In |iV(W 



In 



= c_ #o + lim sup • 



Jo 1 f(u)e- z - r 1,3e ^du 



The integral appearing in the numerator has a leading behaviour given by Lemma 1 



(69) 



f(u)e 



du 



r(q + l) 

r (a+l)/3 



OO, 



the hypothesis (0 — 7r)/3 g] — 7r/2,7r/2[ being indeed satisfied. This shows that the last 
term in the right hand side of (69) vanishes; therefore we have obtained: 



hU\ v Wire*)] a 
h[<p) = limsup — - — = 6*o cos 



7T 



0e 



7T 37T 

2'T 



(70) 



Let us now examine the region e]0, |[. This time — Co < but c + > remains still 
bounded by c_. So using again the relation 6{u) = — 6{u) in the term involving a_ and 
a + we can write this time 



N(z) = e z - rl/l% \J(r,9) 



with 



J(r, 6) = a_ Jq 1 f(u)e~ z - rl/Se ^ du + a + e^ z — z +> 1/3e " f(u)e~ z + rl/H ^ du 

+a e- z - rl/3e " Jo f(u)e- zorl/3 °W du . 

The absolute value of the last two integrals vanishes exponentially for large r so we have 
again the relations (69), (70) from which we conclude similarly that 



h(4>) = 9 cos 



— 7T 



o,- 

' 2 



The analysis for 



[37T 



,2ir[ is similar to the one of the previous interval, up to 



an interchange of — c and c + . So we have obtained the Phragmen-Lindelof indicator for 
e]0, 2n[. We then extend our result to [0, 2ir] using its continuity property and this ends 
the proof. □ 

Proposition 13 proves that all the Nevanlinna matrix elements have the same indicator, 
which we state as 



Proposition 14 The Nevanlinna matrices given by Propositions 8 and 10 have one and 
the same Phragmen-Lindeldf indicator 

= O cos , 0e[O,24 

It follows from relation (57) that they have one and the same type a = 9 . 
Remarks : 

1. Our results agree with the general theorems proved in [2] according to which, for a 
given Nevanlinna matrix, all of its four elements have the same order, type and Phragmen- 
Lindeldf indicator. These quantities are therefore intrinsic to a given matrix. Our Proposi- 
tions 11 and 13 allow us to deal, in one stroke, with the two different Nevanlinna matrices 
obtained. 

2. It may be interesting to observe that the quantities characterizing the asymptotics 
in the complex plane of the Nevanlinna matrix are independent of the parameter c. 

3. It was shown in [4], that for c = 0, the leading term in the asymptotic behaviour 
of the N-extremal mass points was given by 

Since the order is always 1/3 no matter what c is, the leading cubic dependence on n will 
remain, but there is still an interesting open issue: will the coefficient in front of n 3 remain 
independent of c ? 
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